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We study halo nuclei using a two-frequency shell-model approach employing wave functions of 
two different oscillator constants huj in and hui out , the former for the inner orbits and the latter for the 
halo (outer) orbits. An initial application has been made for the halo nuclei 6 He and 6 Li, with Osi/2 
taken as the inner and ( Ops/2, Opx/2) as the halo orbits. Starting from the Paris NN interaction, 
we have derived a G-matrix folded-diagram effective interaction for this two-frequency model space, 
using an essentially exact treatment for the Pauli exclusion operator for the G-matrix. While keeping 
huiin fixed, we have performed calculations with different choices for Tuo ut, treating it as a variation 
parameter. For hiv in = 19.7MeV and hu out — 8.2MeV, our calculated valence energies for 6 He and 
6 Li are -2.77 and -3.55 Mev, respectively, both in good agreement with experiments. The importance 
of certain three-body-force diagrams is discussed. 
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I. INTRODUCTION 



Since several years ago, radioactive-beam nuclear physics has been rapidly developing and has attracted the attention 
of many nuclear physicists With radioactive-beam accelerators, one is now able to measure the properties of 

a large number of drip-line nuclei, which are far from the nuclear stability valley. These nuclei are generally barely 
bound (or unbound) and their structure is typically that of a tightly bound inner core with a few outer nucleons which 
are loosely attached to the core. For example, the halo nucleus 6 He is presumably made of a 4 He core with a pair of 
loosely-attached outer nucleons. The outer nucleons are spatially extended and have a large rms radius, forming a 
"halo" cloud rimming the nucleus as illustrated in Fig. 1. Such nuclei are usually called "halo nuclei". Vaagan et al. 
[PJ have suggested a more distinctive name for them, namely "Borromean nuclei", as a number of halo nuclei, such as 

Li, are believed to form a special type of bound three-cluster system. 

Studies of halo nuclei, or nuclei far off stability, are of interest in many ways. In astrophysics, the location of the 
r-process (rapid proton and/or neutron absorption process) path plays a crucial role for nucleosynthesis. The r-process 
consists of a chain of nuclear reactions, involving a number of intermediate nuclei which are far off stability. To locate 
the path of the r-process accurately, it is necessary to explore, as much as we can, the various properties of such nuclei. 
With radioactive beam facilities, the cross section for certain such nuclear reactions important to nucleosynthesis can 
now be directly measured. For example, the cross section for 8 Li(a, n) 11 B has been measured at the radioactive beam 
facility of RIKEN g. 

Halo nuclei may play also an important role for the synthesis of super-heavy elements, which has been pursued 
by many physicists. The use of halo nuclei, which are already far off stability, as projectiles may provide a new and 
possibly more efficient avenue for producing super-heavy elements. Also from a theoretical viewpoint, the study of 
halo nuclei is of high interest. The empirical nuclear mass formula is, so far, based on a fit to the observed properties 
of stable nuclei. We need information about halo nuclei to generalize the mass formula to regions far off nuclear 
stability. 

It seems to be still rather uncertain about what would be an appropriate and convenient theoretical framework for 
the description of halo nuclei. As will be discussed in Section 2, the ordinary nuclear shell model may not be suitable 
for halo nuclei. We would like to suggest a two-frequency shell model, using single particle (sp) wave functions of 
two different oscillator constants %u>i n and h,Lo out , the former for the inner orbits while the latter for the halo (outer) 
orbits. 

It is our great pleasure and honor to dedicate this short article to Professor G.E. Brown, or Gerry as known to all 
his friends, colleagues, students and tennis partners, on the occasion of his 70th birthday. Gerry has pioneered the 
G-matrix approach for microscopically deriving the nucleon-nucleon effective interaction (V e //) in nuclei. O] 

This approach has turned out to be quite successful, for ordinary nuclei. For example, the CalTech group has 
recently carried out shell-model Monte-Carlo calculations |l3j using G-matrix effective interactions, with results in 
remarkably good agreement with experiments over a wide range of nuclei. In this article, we would like to discuss an 
extension of the above G-matrix approach to halo nuclei, and report some preliminary results. 

In the following, we shall first discuss our two- frequency model space and G-matrix together with some calculational 
details in Section 2. The importance of certain three-body diagrams to the effective interaction will be addressed. In 
Section 3, we shall present some preliminary results for halo nuclei 6 He and 6 Li. A summary and discussion will be 
presented in Section 4. 



II. TWO-FREQUENCY MODEL SPACE AND G-MATRIX 

The rms radius for a nucleon in a harmonic-oscillator shell-model orbit ip n i is given as 

r rms (nl) = (2n + l + ^) 1/2 &; with b = (Ti/muj) 1/2 . (1) 
To fit the observed nuclear radii, there is an empirical formula for choosing the oscillator parameter hu>, namely 

hu = 45A- 1 / 3 - 25A- 2 / 3 (2) 

where A is the nuclear mass number. Using this formula, the empirical value of tuv is approximately 17.2 MeV for 
6 He, corresponding to a length parameter b= 1.55 fm. 

In ordinary shell model calculations, it is customary to use the above empirical formula to determine the fuv value 
used for the shell- model wave functions. For example, we may use hu>=14 MeV (b= 1.72 fm) for most of the sd-shell 
nuclei. The situation for halo nuclei is, however, different. For instance, the empirical rms radius for 6 He is about 2.5 
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fm [Q, which corresponds to a b value much larger than the above value of 1.55 fm. For halo nuclei such as 6 He, we 
need to have wave functions with larger spatial distribution, compared with the ordinary shell-model wave functions. 

Let us also look at the binding energies of some nuclei in the vicinity of 6 He. The observed values JlH] are: 28.29 
( 4 He), 27.40 ( 5 He), 26.32 ( 5 Li), and 29.26 ( 6 He), in MeV. We see that the two-neutron separation energy for 6 He is 
only about 1 MeV, and the single-nucleon separation energies for 5 He and 5 Li are both negative (unbound). It is 
clearly seen that, comparing to ordinary nuclei, 6 He is a very "special" breed. It appears to have a tightly-bound 4 He 
core with a pair of very loosely bound outer nucleons. For this situation, it would be very difficult for the ordinary 
shell model, where all the wave functions have a common hu value, to give an adequate description for the nuclear 
wave function; to have a good description for the 4 He core we need to use shell-model wave functions with a small b 
value, ~ 1.45/m, which is much smaller than the b value needed for a good description for the outer nucleons. In fact, 
to avoid this serious difficulty, many physicists (see [1] and references quoted therein) have abandoned the ordinary 
shell-model approach, and employed instead the method of hyperspherical harmonics jjj for the description of halo 
nuclei. 

We would like to consider a two-frequency shell-model (TFSM) approach for halo nuclei, as described below. In 
recent years, advance has been made for performing realistic nuclear structure calculations starting from a fundamental 
nucleon-nucleon (NN) potential, such as the Paris potential |l6|]. The effects of correlations, which are not explicitly 
included in the model space, are taken into account by deriving the effective hamiltonian for this model space from 
the realistic NN interaction. For the hyperspherical harmonics method, there is still the restriction that one has to 
employ a phcnomcnological effective interaction, such as a Gaussian interaction. This seems to be a drawback. To our 
knowledge, it is not yet clear how to incorporate realistic NN interactions into the hyperspherical harmonics method. 
In contrast, for our TFSM approach it seems to be straightforward to perform microscopic calculations for halo nuclei, 
starting from realistic NN interactions. Otsuka and his coworkers have proposed a variational shell- model approach 
for treating halo nuclei, and have carried out extensive calculations using the Skyrme effective interactions VSM and 
SIII. 

Let us use the p-shell nuclei 6 He and 6 Li to illustrate our approach. It is reasonable to consider the "core" of these 
nuclei as an ordinary alpha particle. Hence we take in the present work the Osi/2 orbit as an inner orbit with oscillator 
constant hui n — 19.7 MeV, corresponding a length parameter 6^ = 1.45 fm. The outer, or halo, nucleons of these nuclei 
are spatially extended, and for them we use the halo orbits Op^/2 and 0pi/2 with oscillator constant hu) ou t, which is 
treated as a variation parameter. In this way, our sp model space may be written as 

P sp ={cj) in ,d in ; cj) out ,d out } (3) 

with the total dimension of P sp being d = (d rn +d out ). The superscripts in and out refer to the inner and outer orbits, 
respectively. Clearly one needs to have all the sp wave functions {4> ln , <p out } be orthonormal to each other, and one 
may be concerned as to how can one fulfill this condition, as htUi n is in general not equal to hio ou t- For a small model 
space, the orthogonality condition is easily satisfied because of parity and angular momentum conservation. For a 
general case, we need to use a common huj value for all the orbits with the same I and j values. Thus, for our present 
6 He and 6 Li case, we use hum not only for the Osi/2 orbit but also for all the other s\/2 orbits, in order to fulfill the 
above orthonormality condition. (Note that the ns\/2 orbits, with n > 0, are needed for the intermediate states in 
the evaluation of the effective interaction as discussed below.) 

For a general model space, the model-space effective interaction V e f / can be derived formally from a folded-diagram 
method |17| ]. Usually this derivation is carried out in the following three steps. First, we calculate the model-space 
Brueckner G-matrix defined by the integral equation fl8| , |l9|| 

G{u>)=V + VQ 2 -L_Q 2 GM, (4) 

w - Q21 Q2 

where to is an energy variable. Q2 is a two-body Pauli exclusion operator, and its treatment is very important as 
we shall later discuss. T is the two-nucleon kinetic energy. Note that our G-matrix has orthogonalized plane-wave 
intermediate states. In the second step, the irreducible vertex function Q — box is calculated from the above G-matrix. 
Finally the energy-independent effective interaction is given by the folded-diagram series 



V eff = Q-Q J Q + Q JQJQ-Q JQJQJQ- -, (5) 

where J denotes a generalized fold, and Q and Q represent the Q — box ]l7| . 

The above formalism is general, and is applicable to either the situation with a single oscillator frequency as in 
earlier shell-model calculations, or the present situation with two oscillator frequencies. For the latter, the calculation 



of the G-matrix is more complicated. The exact solution |18| of the G-matrix of Eq.(4) is 
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G = Gf + AG (6) 

where the "free" G- matrix is 

G F (u) = V + V-^-=G F {w), (7) 



U3-T 



and the Pauli correction term AG is given by 



AGH = -G f (uj)-P 2 , . 1 P 2 -G f H (8) 

e P 2 [l/e + (l/e)G F (w)(l/e)JP 2 e 

where e = lu — T. The projection operator P 2 is defined as (1-Q 2 ). 

The basic ingredient for calculating the above G-matrix is the matrix elements of G F , the free G- matrix, within 
the P 2 space. This space is composed of, however, wave functions of two frequencies, hu>i n and hu> out . This poses a 
technical difficulty because transformations from two-particle states in the cm. coordinates to those in the laboratory 
coordinates are not as easy to perform as in the case of one oscillator frequency. We have adopted an expansion 
procedure to surmount this difficulty, namely expanding the oscillator wave functions with fuuJi n in terms of those 
with tiuJout, or vice versa. When huji n and huj ou t are close to each other, this procedure is relatively easy to carry out. 
But when they are significantly different, the two- frequency G-matrix is considerably more complicated to calculate 
than the one-frequency one. 



We write the projection operator Q 2 as 



Q2 = Q( ab )\ab)(ab\, (9) 



all ab 



and define Q(ab) as 

Q(ab) = < 



0, if b < m, a < 773; 
0, if b < n 2 ,a < n 2 ; 
0, if b < 113, a < ni; 
, 1, otherwise. 

As shown in Fig. 2a, the boundary of Q(ab) is specified by the orbital numbers (711,712,713). We denote the shell 
model orbits by numerals, starting from the bottom of the oscillator well, 1 for orbit Osi/2, 2 for 0p 3 / 2 ,...7 for 0/7/2 
and so on. m denotes the highest orbit of the closed core (Fermi sea). n 2 denotes the highest orbit of the chosen 
model space. We consider here 6 He and 6 Li with 4 He treated as a closed core, thus we have n\=l. Suppose we use a 
model space including all the 6 orbits in the s, p and sd shells. Then for this case n 2 —6. In principle one should take 
713 = 00 Jl8| . In practice, this is not feasible, and one can only use a large 713 determined by an empirical procedure. 
Namely, we perform calculations with increasing values for 773 until numerical results become stable. As illustrated 
in Table 1, we see that there is a convergence behavior and for our present calculation a choice of 773=21 appears to 
be adequate. (This 77,3 value will be used in our present work.) 

There is an important point about whether we use a large 773 or not. As illustrated in the diagram of Fig. 3 (upper 
part), we used railed lines to denote nucleons outside the model space while for those inside we use bare lines. This 
diagram has one line outside the model space, and hence is a legitimate G-matrix diagram, even if the state b refers 
to a single-particle state below the Fermi energy (i.e. 6 < n\). We are, however, excluding this diagram from the 
G-matrix when we use a large 773. How can we do this? 

There are 3-body components in the effective interaction. Consider the example of the diagram of Fig. 3a (lower 
part), which represents a 3-body- force (irreducible) diagram, where j denotes, for example, a particle in orbit 0si/ 2 
and a represents an orbit outside the model space. (We now draw diagrams with respect to the bare vacuum.) Note 
that the exchange diagram of (3a) is just diagram (3b) with a minus sign. Diagram (3b) is, however, equivalent to the 
diagram of Fig. 3 (upper part). Hence the diagrams of the type shown in Fig. 3 (upper part) are exactly canceled by 
the corresponding 3-body- force diagrams of the type shown in Fig. 3 (lower part). In other words, such 3-body- force 
diagrams can be incorporated into the G-matrix and in so doing they cancel the G-matrix diagrams of the type shown 
by Fig. 3 (upper part). Consequently we should use a large 773. 

It is common to employ a simpler Pauli operator [£0| (and references quoted therein), defined by " Q(ab) = 
1, if (p a + p ) < Nq; = 0, otherwise", with p a = 2n a + l a and similarly for p^. Such a Pauli operator has the shape 
shown in Fig. (2b). Here Nq is typically a small number, such as 2 for the p-shell nuclei. A main advantage of this 
choice is its convenience for numerical calculation. However, this choice for Q corresponds to a small 773, and hence 
the 3-body-force diagrams as shown by Fig. 3 (lower part) are not yet taken into account by the small-773 G-matrix; 
they have to be calculated separately. In contrast, we use a large 773 in our present work, and hence the effect of such 
3-body-force diagrams is already absorbed in our G-matrix. Using a small-773 G-matrix alone may over-estimate the 
nuclear binding energy, as the effect of such 3-body-force diagrams is usually repulsive. 
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III. APPLICATION TO HE AND 6 LI 

Experimental data about 6 He and 6 Li are rather scarce. The valence energies of these nuclei, defined as 

£w( 6 He) = -{BE( 6 Ue) + BE( 4 Ue) - 2 x RE( 5 He)] (11) 

and 

E va i( 6 U) = -[BE( 6 U) + BE( 4 Rc) - BE( 5 Re) - BE( 5 Li)}, (12) 

are nevertheless well known. From the empirical binding energies (BE) [ [TBI , the valence energy for 6 He is found as 
-2.77 MeV and that for 6 Li is -3.83 MeV. 

We want to calculate these quantities using a TFSM approach. For the inner orbit 0s!/ 2 we use &i„=1.45 fm. For 
the outer orbits 0p 3 / 2 and Opi/ 2 we use a sequence of b ou t values, ranging from 1.45 to 2.25 fm. The G-matrix is 
first calculated using a Pauli operator specified by (ni, n 2 , n3) = (l,6,21), starting from the Paris potential. Then, 
following closely the Q — box folded-diagram procedure of Shurpin et al. we derive the matrix elements of the 
effective interaction V e ff for the p-shcll. (In the Q — box we have included diagrams first- and second-order in G. 
This includes the diagram of second order in G with an intermediate state of two nucleons in IsOd shell. To avoid 
overcounting we have to use n 2 = 6 in the Pauli operator for the G-matrix equation.) Our results for the matrix 
elements of V e // are listed in Table 2. As expected, these matrix elements are seen to become generally weaker as b ou t 
increases, a reflection that the mean distance between the interacting valence nucleons is larger. The empirical (6-16) 
matrix elements of Cohen and Kurath are also listed for comparison. For some cases, the agreement between our 
calcualated matrix elements and their values becomes improved when we use a 6 out value larger than 1.45 fm. 

To obtain the valence energies, we just diagonalize the two-particle matrix (Hq + V e ff) within the p-shell, using 
the matrix elements of Table 2. Hq is the sp Hamiltonian. For the valence-energy calculation, the single-particle 
energy for p 3 / 2 is defined as zero. But that for pi/% is uncertain, both experimentally and theoretically. For an initial 
calculation, we have fixed it as 10 MeV. (Using an empirical formula given in Bohr and Mottelson's book p2| , the 
P1/2 — Pi/3 spin-orbit splitting is estimated to be ~ 10.2MeV.) The valence energy for 6 He is then just the lowest 
eigenvalue of the T=l, J=0 secular matrix, and that of the T=0, J=l matrix gives the valence energy for 6 Li. Our 
results for these energies are displayed in the lower part of Fig. 4, calculated as a function of b ou t- At b ou t — 1.45 
fm, the calculated values are too large. At b out — 2.25 fm, our results are -2.77 and -3.55 MeV for 6 He and 6 Li, 
respectively, both in reasonable agreement with the corresponding empirical values -2.76 and -3.83 MeV. 

We have also calculated the single-particle separation energy of 5 He, defined as E sp — BE( 5 He) — BE^He). Our 
calculation was done using a one-body Q — box consisting of the three 1-body diagrams (first- and second-order in 
G) of Rcf. |H| and then summing up the folded-diagram series consisting of the one-body Q — box. These results 
are presented in the upper part of Fig. 4. There is an improvement of our results as b ou t increases. However, our 
calculated separation energy is still too large at b out ~2.25 fm. It is of interest to note that there seems to be a 
saturation point for E sp at b ou t ~2.25 fm. We recall that at this b out our results for E va i are also in good agreement 
with experiments. Assuming a pure two-frequency s 4 p 2 wave function with bi n =lA5 and b ou t— 2.25 fm , the rms 
radius for 6 He is obtained, using Eq.(l), as 2.37 fm, which is also in reasonable agreement with the empirical value of 

Encouraged by the above preliminary results, we have begun to tackle a more difficult problem, to calculate the 
total binding energies of He, 5 He and 6 He using our present two-frequency approach. Here we employ a two-shell ( 
s and p) model space, treating all nucleons of the nucleus as active. For s-shell we use fixed length parameter 6j„ = 
1.45 fm. For p-shell we use a varied b out . The components of a spurious cm. motion, which are contained in such 
a model space including valence orbits of 2 major shells are removed using the technique of Gloeckner and Lawson 
j23j . Our preliminary results indicate that for 4 He an energy minimum seems to be at b ou t — 1.45/m, while for 
°Hc the minimum seems to be located at 6 ont ~ 2.25/m. Our calculated total binding energies have turned out to 
be too small, for both nuclei, compared with experiments. This is perhaps a common problem of nuclear structure 
calculations with no-core employing realistic NN interactions. Further studies are being pursued, and we hope to 
report our results before long. 

IV. SUMMARY AND DISCUSSION 

Halo nuclei have a special nuclear structure; they have typically a tightly bound core with a few halo nucleons 
loosely attached to the core. One would encounter a lot of difficulties, in describing such nuclei using the conventional 
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shell model where the basis functions all have the same hto value. For example, the wave functions with 6j n =1.45 
fm would give a good description for the alpha-particle core of 6 He, but they are spatially much too compressed to 
describe the halo neutrons of 6 He. If one insists to use the compressed wave functions also for the halo nucleons, 
then one needs to use a very large active space (such as a space with s, p, sd and pf shells for the case of 6 Hc). This 
would be very tedious as well as not economical. A much simpler solution is to allow for a two-frequency model space, 
where one uses 6j„ for the core and b ou t for the halo part of the nucleus as described in the present work. We have 
obtained some preliminary results using this approach. Our results seem to indicate that this two-frequency method 
is a feasible and promising approach for halo nuclei. 

What effective interactions should one use in a two-frequency shell-model calculation? It would be very difficult 
to determine them by an empirical best-fit procedure, which has been successfully used for one-frequency shell-model 
calculations p5| . This is because now the length parameter b out is a variable and the matrix elements of the effective 
interaction are dependent on it. A basic and rather "ambitious" solution to this problem would be to derive the 
effective interaction directly from a fundamental NN interaction. We have carried out some preliminary studies in 
this direction. We first calculate the two-frequency G-matrix, with an accurate treatment for its two-frequency Pauli 
exclusion operator. Then the effective interaction is obtained via a standard folded-diagram method. Our calculated 
valence energies for 6 He and 6 Li are in satisfactory agreement with experiments. As indicated by our preliminary 
results, 4 He seems to favor a one-frequency (bi n =b out =lA5 fm) shell model, while a two-frequency (6j„=1.45, b ou t=2.25 
fm) one is perhaps needed for 6 He, 6 Li and 5 He. 
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-20 
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3 


-16.048 


-15.424 


-14.453 


11110 1 


6 


-15.295 


-14.706 


-13.842 


11110 1 


15 


-15.238 


-14.637 


-13.767 


11110 1 


21 


-15.238 


-14.637 


-13.767 


2 3 2 3 2 


3 


-6.138 


-5.691 


-5.132 


2 3 2 3 2 


6 


-5.635 


-5.339 


-4.921 


2 3 2 3 2 


15 


-5.563 


-5.277 


-4.877 


2 3 2 3 2 


21 


-5.563 


-5.277 


-4.877 


111110 


3 


-9.456 


-9.360 


-9.213 


111110 


6 


-9.151 


-9.093 


-8.995 


111110 


15 


-9.140 


-9.083 


-8.986 


111110 


21 


-9.140 


-9.083 


-8.986 



Tabic 1. Dependence of the two-frequency G- matrix on the choice of n^. Listed are the matrix ele- 
ment (abTJ\G(oj)\cdTJ), in MeV, calculated with the Paris potential. The orbits 1,2,3 represent, respectively, 
0si/2, 0p 3 / 2 , 0pi/2- We have used bi n ~ 1.45 and b out — 2.0 fm for the length parameters, and m = 1 and n 2 = 3 
for the exclusion operator. 
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1.45 


1.75 


2.00 


2.25 




2 2 2 2 1 


-1.76 


-2.15 


-2.28 


-2.24 


-3.14 


2 2 2 3 1 


5.64 


4.69 


3.89 


3.18 


4.02 


2 2 3 3 1 


2.55 


2.24 


2.13 


2.03 


1.09 


2 3 2 2 1 


5.81 


4.64 


3.78 


3.05 


4.02 


2 3 2 3 1 


-7.63 


-6.75 


-5.96 


-5.15 


-6.54 


2 3 3 3 1 


3.12 


1.72 


0.83 


0.20 


1.39 


3 3 2 2 1 


2.72 


2.20 


2.01 


1.87 


1.09 


3 3 2 3 1 


3.11 


1.82 


1.08 


0.56 


1.39 


3 3 3 3 1 


-3.13 


-2.80 


-2.61 


-2.37 


-4.26 


2 2 2 2 1 


-3.43 


-3.28 


-2.92 


-2.51 


-2.74 


2 2 3 3 1 


-5.21 


-3.88 


-3.02 


-2.38 


-5.32 


3 3 2 2 1 


-5.41 


-3.86 


-2.99 


-2.36 


-5.32 


3 3 3 3 1 


0.44 


-0.45 


-0.74 


-0.81 


0.34 



Table 2. Some p-shell effective interactions calculated with various b ou t values. Listed are the matrix element 
(abTJ\V e ff\cdTJ), calculated from the Paris potential. The orbits 2,3 represent, respectively, 0p 3 / 2 ,0pi/ 2 - We have 
used bi n — 1-45 fm and bare G-matrix calculated with (ni,n2,n^)— (1,6,21). The empirical (6-16) matrix elements of 
Cohen and Kurath |^] are listed under column CK. 
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FIG. 1. A halo nucleus. 




FIG. 2. Pauli exclusion operator Q2- 
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FIG. 3. A second-order diagram belonging to the G-matrix (upper part) and Three-body- force diagrams (lower part). 
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FIG. 4. The lower part of this fi gure shows the valence energy E va i for He and Li as a function of b out , while the upper 
part displays the single-particle separation energy E sp for He as a function of bout- 
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